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O . Abstract 
<N 

\& \ A new multiplicity distribution with multifractal properties which can be used 

in high-energy physics and quantum optics is proposed. It may be considered as 
. a generalization of the negative-binomial distribution. We find the structure of 

the generating function for such distribution and discuss its properties. 

^ 1. Introduction 

^3 ■ 

Multifractal analysis in high-energy physics and quantum optics has received a great 
interest in recent years due to the possibility to obtain quantitative and qualitative 
results concerning multiparticle production in different processes. In addition, the 
analysis is becoming an important theoretical tool to discriminate between different 
dynamical models. 

The fact that the normalized factorial moments (NFMs) 

(w(w-l)...(w-g+l)) 
{n) q 

(n is the number of particles in a restricted phase-space interval 5, (. . .) is the aver- 
age over all events) depend on the size of the phase-space interval (bin) as F q ~ 5~^ q 
(intermittency phenomenon) is a manifestation of non-statistical fluctuations in the 
multiplicity distribution of secondary particles produced in high-energy physics [|T]-0|. 
The multifractal behavior, when the anomalous fractal dimension d q = <fi q /{q — 1) 
depends on q, is particularly important because such behavior is typical for the vast 
majority of experiments [EH4J]- It is more pronounced in two- and three-dimensional 
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phase-space domains. This behavior has also been found in photon-counting experi- 
ments on laser fluctuations near threshold, where S is the counting time interval T || . 
Thus, the problem of multifractality is a common one both for high-energy physics and 
quantum optics. 

In this paper, we shall discuss the theoretical aspect of the problem in the context 
of high-energy physics by means of bunching parameters (BPs) |J. We introduced 
this quantity in order to get a simple and efficient method for the analysis of complex 
multiplicity distributions in restricted phase-space intervals. Our consideration can 
also be used in any field of research, where local dynamical fluctuations are a subject 
of investigation. 

As is well known, the negative-binomial distribution has become the focus of interest 
in view of its applicability to the study of multiparticle productions in high-energy 
physics. However, it has been noted that this distribution has no multifractal properties 
for small phase-space intervals 

Moreover, the study of charged-particle multiplicity distributions in restricted ra- 
pidity intervals conforms that the negative-binomial distribution is not sufficient to 
describe the data in Z° hadronic decays due to a shoulder structure of the experimen- 
tal distributions ||. A similar conclusion was obtained by ALEPH collaboration ||: 
it was shown that the negative-binomial distribution does not describe the experimen- 
tal data, either in restricted rapidity intervals or in the full phase space. The UA5 
collaboration has also observed that the negative-binomial distribution fails to give a 
good fit to data at the center-of-mass energy of 900 GeV |T(| in full phase space due to 
a shoulder structure. This structure is explained by the superposition of 2-jet events 
with low multiplicities and 3-jet and 4-jet events yielding much larger multiplicities. 
As a rule, this leads to more complex multiplicity distribution for full phase space than 
the negative-binomial one. 

The conclusion must be that the negative-binomial distribution is not sufficient 
to describe the experimental distributions both for restricted rapidity windows and 
for full phase space in definite experimental situations and, hence, a true multiplicity 
distribution must be more complicated. 

In this paper we propose a new multiplicity distribution which has multifractal 
properties for small phase-space intervals and is very similar to conventional negative- 
binomial form for large phase-space intervals. We shall analyze this distribution in 
terms of BPs and bunching moments. As we shall see, such an investigation is sim- 
pler than the analysis of multiplicity distributions with the help of NFMs. Moreover, 
recurrence relations for probabilities can lead to a non-traditional form of generating 
functions, both for well-known distributions (Poisson, geometric, logarithmic, positive- 
binomial, negative-binomial distributions) and the multiplicity distribution with mul- 
tifractal properties for the small phase-space intervals. 

In the second section, we give a short introduction to the bunching-parameter 
method. In the third section, we consider the general form of a multifractal negative- 
binomial like distribution and derive its generating function. In Sect. 4, the properties 
of Markov branching process leading to this multiplicity distribution are considered. In 
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Sect. 5, we discuss a particular form of such a generalization of the negative-binomial 
distribution and illustrate its multifractal properties. In Sect. 6, we present the con- 
clusions. 

2. Bunching-parameter approach 

Normalized factorial moments have become an important and popular topic of 
experimental and theoretical investigations in high-energy physics and quantum optics. 
Measuring the NFMs is equivalent to measuring the multiplicity distribution. On the 
other hand, recently another simple mathematical tool to investigate the behavior of 
the multiplicity distribution in different phase-space intervals has been proposed. In 
order to reveal "spike" structure of events, one can study the behavior of probabilities 
themselves by means of the BPs. The definition of the BPs is given by the formula |J 

= ~ 1 p2 > Q > 2, (2) 

where P q is the probability of finding q particles inside the limited phase-space interval 
5 0. For example, T] g — 1 for the Poisson distribution. If the size of the phase space is 
small and the average number of particles in this interval is approximately proportional 
to 5, then we have the following approximate relation between the NFMs and the 
BPs § 

F^flvr 1 " 1 - (3) 

Such relation, in fact, means that BPs share with the NFMs the property of suppression 
of Poissonian noise. 

For the intermittent fluctuations, one expects 772 oc 5~ d ' 2 (c^ is the second-order 
anomalous fractal dimension), while the high-order BPs may have a different depen- 
dence on 5. In it has been shown that for the high-energy collisions with the mul- 
tifractal behavior of the NFMs d q = d 2 {l — r) + d,2rq/2, the power-like behavior of the 
BPs r]2 oc 5~ d2 , 7] s = 7/2, s > 2 is valid. The positive parameter r can be interpreted as a 
degree of multifractality (for r = 0, we have exactly monofractal behavior). Thus, the 
problem of the multifractal multiplicity distribution with inverse power 5-dependence 
of all BPs is central one. 

As discussed in 0, the use of the BPs can give a general answer to the problem 
of finding a multiplicity distribution leading to observed intermittency. Indeed, any 
multiplicity distribution can be expressed as 

On the contrary, a multiplicity distribution can not be expressed in terms of its factorial 
moments if this distribution is not truncated one HI. 
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3. Recurrence relation for multifractal multiplicity distribution 

A necessary and sufficient condition for a multiplicity distribution P n to be negative- 
binomial one (Pn) is the recurrence relation 



^L =tVl = £±^Zi), „>,, (5) 

r n-l ' i 

where a, 6 are positive constants for fixed 5 (6 < 1). Note that throughout this paper 
we shall treat the probabilities and the parameters C n _i, a, 6 as continuous functions 
of the phase-space interval of size 5. If a ^ 0, 6 ^ 0, iterations of (Q) with the 
normalization condition Y^n°=o Pn = 1 & ve the negative-binomial distribution. In the 
limit b 0, the recurrence relation gives a Poisson distribution. The case a — > 
at constant 6 shows that the negative-binomial distribution reduces to a logarithmic 
distribution. For a = b, we get a geometric multiplicity distribution. For a > 0, 6 < 
and a/6 integer, we have a positive-binomial distribution. Using the definition of the 
BPs, (|5]) and theorem of , it is easy to see that none of these distributions leads to 
multifractality if any assumption is chosen for the behavior of a, b in small phase-space 
bins. 

Thus, it is important to find some multifractal generalization of the commonly used 
multiplicity distributions. From a physical point of view, to find such a multiplicity dis- 
tribution, in fact, means to understand the reasons of intermittency with multifractal 
behavior of the anomalous fractal dimensions. However, the level of theoretical under- 
standing of this phenomenon is still insufficient (hadronization problem) and is quite 
different for various types of collision processes |4[] . Nevertheless, from a mathematical 
point of view, we can propose a distribution that has a priori multifractal behavior. 

There is a natural way to include BPs with power-law behavior into a new recurrence 
relation in order to obtain a modification of (^) which can generate the multifractal 
multiplicity distribution in the limit of small 5. To see this, let us rewrite the definition 
(D of BPs as follows 

P q _ q-l P q -i 



•q * 9 ■ (6) 
P q -i 9 q P q -2 1 ' 

As mentioned in Sect. 2, for a multifractal behavior with r = 1, all BPs have the same 
power-law behavior, i.e., in the simplest case, we can write 



r] q = g <x 5 13 . (7) 

The f3 is a positive constant and is taken as a measure of the strength of the multifractal 
effect. After that, by combining and (|5|), we assume the following recurrence relation 

P n o + 6(n-l) n-lP fl -i , , 

-p— = — - — +9 ^rp~' (8) 

where n > 1. Here, in fact, g can be either a new parameter or some combination 
of the parameters a and b (the latter case we shall discuss below). Equation (||) is a 
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sufficient condition to construct a distribution which for a large phase-space interval 
{g <C 1) is very similar to the multiplicity distribution (||) and has the multifractal 
behavior for a small phase-space interval (the value of g is not small). 
The expression (H) can be rewritten as 

P ~ n ~ l a 1 

" CU = ££-(a + &(n- /-!)), (9) 



r n-1 {=Q Tl 



with the sum equal to 



bin - ng)- 1 ^ - g n ) + n + (g n - 1)(1 - g)- 1 ], for g ^ 1, 
6(jfe + 0.5(n-l)), for = 1, 



(10) 



where k~ x = b/a was called aggregation coefficient by Giovannini and Van Hove [|TTJ for 
the usual negative-binomial distribution. For a compact description of the recurrence 
relation, we shall use expression for C* n _i in the form of sum (Q). Using @ and @, 
we obtain for BPs 

1 + k-\q- 1) 

TU9 l {l + k-i{q-l-2)) 

We have multifractal behavior because the second term in ([TI]) has power-law behav- 
ior and, for enough small 5, can be the leading one. The multiplicity distribution 
corresponding to (|^), ([]) is 

p n—l s 

Pn = ^l[T,9 l (a + b(s-l)). (12) 

n - s=0 1=0 

Throughout this paper, we consider the probability Po = 1 — Yl™=i Pn of having no 
particles in 5 as a normalization constant. Since in the limit g — > expression (|l2|) 
reduces to a negative-binomial distribution, we shall call this distribution a multifractal 
negative-binomial distribution (MNBD) and denote it as Pff . 

Let us note that we can analyze a multiplicity distribution written in terms of the 
recurrence relations using the generating function for P n /P n -\. Let 

oo ~n p 

""^M^ (13) 

be the generating function for the ratios P n jP n -\. Then, the BPs are given by 



^ = (14) 

For example, the generating function of the negative-binomial distribution has the 
following form 

G N (z) = (a-b)(e z - 1) + bze z . (15) 
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For the Poisson distribution, we have 

G p (z) = a{e z - 1). (16) 

In terms of the generating function, the recurrence relation (^j can be rewritten in 
form of the following differential equation 

^M = ^l + g G^, (17) 

dz dz 
with the initial condition (see (|T3|)) 

G M (z = 0) = G N (z = 0) = 0. (18) 

Using (|17D, we can obtain the generating function in integral form as follows 

CO „ 

G M (z) =G N (z) I G N {z)d i z. (19) 

i=i 

Using ( |15D and condition (|l^), one gets the analytical solution of the equation (|17|) 

"«-I^(' + »-T^)-' ! ^ i - (20) 

for g 7^ 1 and, using (|9|) and (^3|), G M (z) = e z (az + bz 2 /2) for g — 1. Below, we shall 
not consider the trivial case, when g — 1. 

Now let us mention two limiting cases: 

a) In the simplest Poisson case, when b = 0, we have the following multifractal 
Poisson distribution with the generating function of the form 

G MP (z) = (e z - e z9 ) . (21) 

Then, from (JTT|) and (|T2"D, one has 

n n—1 s 

e P = ^IIE^ (22) 

' t - s=0 Z=0 

<r = =^ + 9. (23) 

Here a is not an average multiplicity as for a usual Poisson distribution. 

b) The limit a — > is also interesting since it leads to the multifractal logarithmic 
distribution 

gMI m^-t^H*(t4f (24) 
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From ([TI]) and fll2|), one gets 



Un—l n—2 s 

Pn L = Pi-rXlY,9 l {s-l + l\ (25) 



u - s=0 1=0 



< L = ^-2_",_ , ^ + 9, (26) 



where q > 2 and Pi is a normalization constant (for the logarithmic multiplicity dis- 
tribution, Pq = 0). We see that these multiplicity distributions have the same power 
law behavior of the high-order BPs for small 5. In this sense, the distributions are 
equivalent for small phase-space intervals. 

It is important to emphasize here that an infinite sequence of probabilities P n can 
be normalized if it converges, i.e., if P n — > for n — > oo. This is possible if, for every i 
greater than some number £, we have the following condition 

P ' <1 (27) 



Pi 

(ratio test). For the MNBD © this is possible if 

0<o<l, 0< < 1. 28 

1-9 

For any other domain of g, we must truncate the MNBD, putting Pf 1 = for i > 

One can understand, from definition (0), that BPs are sensitive to the local multi- 
plicity fluctuations (or to the behavior of the multiplicity distribution in small phase- 
space intervals) near the multiplicity q — n—1. In order to study the total contribution 
from multiplicity fluctuations for large values of n, it is appropriate to introduce the 
"bunching moments" b q as follows 

^Ml— g(7^b < 29 > 

following an analogy with factorial moments. The knowledge of the generating function 
gives us a possibility to calculate both the BPs and the bunching moments. The 
higher the rank of the b q , the more sensitive they are to the "tail" of the multiplicity 
distributions for large n. In order to normalize the bunching moments B q , we put 

b I i 

B q - f = 7,rM (30) 

Oi Cr | z= i 

Then for the Poisson distribution (0) we have 

5f = l, g=l...oo. (31) 
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For the negative-binomial distribution (|T5|): 

B? = $±4, g=l...oo. (32) 



This means that the negative-binomial distribution is broader than the Poisson dis- 
tribution (B^ > Bq). For the geometric multiplicity distribution (k = 1), the nor- 
malized bunching moments are larger than those of the negative-binomial distribution 
with k > 1. In this case, we can say that the geometric distribution is broader than 
the negative-binomial one. For the positive-binomial multiplicity distribution, where 
k < and integer, the normalized bunching moments are smaller than unity because 
this distribution is narrower than Poisson. 
For MNBD fl20[), one gets 

R M = K9 - i)^ 9 ' 1 - i) + g q e g - 1 -g(g + i) + g r „v 

k(g-l)(ge^-l)+ge^-2g + l " ,0 °' 



q = 1 . . . oo. For a small g (or for a large phase-space interval), the MNBD slightly 
differs from the negative-binomial distribution. For g — > (B^ 1 — > -B^), the MNBD 
tends to the negative-binomial distribution. The situation drastically changes when 
large g (or small 5) is considered. The B^ 1 becomes larger than B^ , since the MNBD 
becomes broader than the negative-binomial distribution for all k. Bear in mind of this 
property, it is not surprising that for this region of phase space the MNBD reflects the 
increase of the local intermittent fluctuations and can lead to a multifractal behavior 
of the anomalous fractal dimension. 

The generating function Q = J2™=o z n P n commonly considered in probability theory 
and its applications can also be used in the multifractal generalization ||12|| . To analyze 
this generating function, it is appropriate to use the usual NFMs or the normalized 
cumulant moments [[|. In fact, these two "languages", using the BPs and the NFMs 
for the study of multiplicity distributions in terms of the generating functions G and 
Q, are equivalent. Nevertheless, in some cases, the analysis of multiplicity distribu- 
tions in small 5 with the help of BPs and the bunching moments is simpler, because 
an analytical form of the NFMs for some multiplicity distributions may be too compli- 
cated. For example, to find a simple form of the generating functions for multifractal 
distributions flT2D , ([22]), (f25|), is rather difficult. The technical advantage of the use of 
these quantities comes mainly from the fact that, as a rule, the structure of the ratio 
P n /P n _i is simpler than the form of the probabilities P n themselves. 

It is appropriate to make some remarks here on the relationship between the form 
of the recurrence relation (^) and the definition of the generating function (|llj). The 
particular form of recurrence relation (|8]) was chosen because of its simplicity. For 
example, the factor (n — l)/n in the last term of (|8|) was chosen only because it is 
possible to rewrite this relation in the form of generating function ( [13]) . It is easy to 
check that other similar forms involving P n -\ and P n -2 in the recurrence relation can 
lead to qualitatively similar results, providing a singular form of BPs for all order. 
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However, in these cases, one needs to introduce other forms of generating function in 
order to obtain a differential equation for the MNBD with simplest solutions. 

4. The MNBD as a solution of non-linear Markov process 

Now we shall show that the form of the MNBD can be obtained from a stochastic 
Markov process with non-linear birth rate. Let P n (t) be the probability to have n 
particles at time t. Of course, such a choice of an evolution parameter is not unique. 
In principle, the evolution variable t can be connected, for example, with the squared 
mass of the branching parton in the parton shower. For simplicity, we shall assume that 
the process starts at time t = 0, with the initial condition P (t = 0) = 1, P n (t = 0) = 0, 
n > 0. We shall consider a very general birth-death process with an infinitesimal birth 
rate w£ and an infinitesimal death rate w~ of particles, treating these parameters as 
continuous functions of t. The corresponding Markov equation is 

P' n {t) = wZ^Pn-i{t) + w- +1 P n+1 (t) - {w+ + w-)P n (t), (34) 

fill . For a stationary process, when time goes to infinity, the P n (t) are t-independent 
constants. Then, from ( |34"D one has n n — 7r n+1 = 0, n — 1, 2, .., ir n — w~P n — Wn-xP n -i- 
Taking into account 7r = 0, we have n n = 0, n > 1 and, hence, a solution of (|34]) in 
the form of the recurrence relation 

Pn w n-l ( r lK\ 

Pn-1 W- ' 1 ' 

The negative-binomial distribution can be considered as a stationary solution of (|34|) 
if we put the linear forms for and w~, 

Wn = 7 + fln, w~ = pn, (36) 

where a = j/p and b = j3/p (7, (3, p are ^-independent ). If we admit that the parameter 
is the non-linear function of n, i.e., 

n 

< = Y,9 l h + P(n-l)}, w~=pn, (37) 
z=o 

then the MNBD can be found as a stationary solution of the evolution equation (|34|). 
Here, in fact, parameter g represents the strength of influence of the non-linearity in 
the equation. For 5 — > the non-linearity of the birth rate increases. 

So, the multifractal structure of MNBD can be caused by the non-linearity of the 
stationary Markov process. From the point of view of high-energy physics, the multi- 
particle production in QCD and the subsequent transition to hadrons have a strongly 
non-linear nature. In this sense, the form of (|37D may reflect the non-linearity of the 
underlying multiparticle dynamics leading to multifractality. 
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5. The choice of parameter g and application 
to experimental data 

Up to now, we have considered g as a free parameter with the power-law behavior 
g oc 5~P. The next question is how does one choose g, if one wants to obtain the 
multiplicity distribution with a degree of multifractality r (see Sect. 2). The simplest 
way is to assume the following form 

/u\ r 

= co(r)k- r , (38) 

where u(r) is some function which tends to zero for r — > (in the simplest case, u)(r) = 
r). The parameter r allows interpolation between the negative-binomial distribution 
with monofractal behavior (r = 0) and the MNBD with the multifractal behavior (r > 
0). Let us remind that the aggregation coefficient k~ l is related to the mean multiplicity 
< n >= a/ (1 — b) and the dispersion D of the negative-binomial distribution as D 2 =< 
n > + < n> 2 

From the point of view of high-energy physics, our choice is justified by the following 
reasons: 

i) If we have the negative-binomial distribution, for a large (pseudo)rapidity interval 
k~ x has a small value (k^ 1 ~ 0.1 — 0.01). Applying a fit by ([12]) , we can hope that g 
in the MNBD will be small also and this distribution will be similar to the negative- 
binomial distribution. 

ii) It is essential that the negative-binomial distribution is approximately valid, not 
only for large phase-space intervals, but also for small ones. For the negative-binomial 
distribution, the behavior k^ 1 oc 5 _/3 lies in the framework of the assumption that 
intermittency is governed only by the aggregation coefficient. In J7| Van Hove showed 
that, if k^ 1 oc 5~" for small S, then the negative-binomial distribution has monofractal 
behavior, F q oc 5^ dqt - n ^ l \ d q = (3. Then r] 2 oc i] s = const., s > 2 p. In our case, 
the assumption u(r) ^ yields intermittency with the multifractal behavior of the 
anomalous fractal dimension. 

From the experimental point of view, a good approximation for the aggregation 



coefficient of the negative-binomial distribution is [14 



r^ctf 2 , M = — , (39) 

o 

where A and 5 are three-dimensional full and limited phase-space intervals, respectively, 
<i 2 is the anomalous fractal dimension of the second order. For different reactions, the 
parameters c and are of the same order of magnitude, moreover, the d 2 has almost 



universal value ~ 0.4 |14|. Expression (p5|) has been obtained from the well-known 
relation = 1 + k~ l between second-order NFM and the aggregation coefficient 
which is correct for the negative-binomial distribution. 

iii) Since the fits for different reactions show a logarithmic increase of k~ x with 
energy ■sfs in the negative-binomial distribution, one may expect that for the MNBD 
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we shall obtain the similar effect. Then, the distribution ([12]) in the full phase space 
may slightly differ from the standard negative-binomial distribution for large energies. 
This is quite important because, as mentioned already in the introduction, the usual 
negative-binomial distribution fails to describe multiplicity distributions at 900GeV in 
pp collisions and in Z° hadronic decay for full phase space. 

Let us obtain the anomalous fractal dimension for the MNBD in the particular case, 
when g = oj{r)k~ r . The calculation of the NFMs can be simplified when an average 
multiplicity in 5 is small. Then, the NFMs are given by expression (^]). Let us discuss 
two domains of the parameter r: 

(i) < r < 1: From flTT|) we have 

7/2 = 1 + k~ x + w(r)k~ r , (40) 

Vs = v - 2 f( LtV+y-u i o\\ + w ^ k ~^ s>2 - ( 41 ) 
22i=o^{r)k rl (l + k L {s-l-2)) 

We see that, if k oc 5 d2 , the leading terms of the BPs have the following behavior: 
r/ 2 oc 5~ d \ T} s oc 5~ rd \ s > 2. Then 

F q oc r^ 1 ) d q = d 2 (l - r) + d 2 r|. (42) 

For r = 0, we have the monofractal behavior d q = d 2 , and the MNBD reduces to the 
negative-binomial distribution. 

(ii) r > 1: The leading terms of the BPs are given by the expression rj q oc 5~ rd2 , 
q > 2. The corresponding anomalous fractal dimension is 

d q = d 2 A. (43) 

The values of r for different reactions have been discussed in . 

6. Conclusion 

We have proposed a new multiplicity distribution with multifractal properties for 
small phase-space intervals, basing on the simplicity of the analysis of multifractality in 
terms of BPs and bunching moments. Guided by the fact that, for the multifractality 
of normalized factorial moments the same power-law behavior for all orders of the 
BPs is necessary and sufficient, in this paper we focus our attention on the analysis 
of the MNBD with the multifractal behavior. The MNBD may be considered as a 
generalization of the negative-binomial distribution with a new free parameter g (or r 
for the particular case, when g = uj{r)k~ r ) which has a power-like behavior for small 
phase-space intervals. 

Theoretically, the question arises what is the physical reason of such multifractal 
distribution (or similar to it). The problem of multifractality is very complicated and 
requires further examination both in high-energy physics and quantum optics. Note, 
as an example, that some versions of cascade models can lead to this distribution 
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because they have the same anomalous fractal dimension (see a-model ||15|| , where the 
anomalous fractal dimension has the form (^) with r = 1). 

On the experimental side of the question, we hope that the MNBD is interesting 
since it yields a new possibility to describe experimental data. 
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